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ABSTRACT 

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation 
models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratios (e.g., B/C) 
produce too few antiprotons. In the present paper we discuss one possibility to overcome these difficulties. Using the 
measured antiproton flux and B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei as 
measured in the heliosphere may contain a fresh local "unprocessed" component at low energies perhaps associated 
with the Local Bubble, thus decreasing the measured secondary to primary nuclei ratio. The independent evidence for 
SN activity in the solar vicinity in the last few Myr supports this idea. The model reproduces antiprotons, B/C ratio, 
and elemental abundances up to Ni (Z < 28). Calculated isotopic distributions of Be and B are in perfect agreement 
with CR data. The abundances of three "radioactive clock" isotopes in CR, 10 Be, 26 Al, 36 C1, are all consistent and 
indicate a halo size z h ~ 4 kpc based on the most accurate data taken by the ACE spacecraft. 



INTRODUCTION 

The result of our recent analysis (iMoskalenko et al. , 200 lbl 2002 ), in agreement with calculations of other authors 
(hvlolnar and Simonl 2001 ), was that matching the secondary/primary nuclei ratio B/C using reacceleration models 
leads to values of the spatial diffusion coefficient apparently too large to produce the required p flux, when the prop- 
agated nucleon spectra are tuned to match the local p and He flux measurement s. Assuming the curren t heliospheric 
modulation models are approximately right, we have the following possibility (IMoskalenko et all 12003 ) to reconcile 
the B/C ratio with the measured flux of secondary p's. The spectra of primary nuclei as measured in the heliosphere 
may contain a fresh local "unprocessed" component at low energies thus decreasing the measured secondary to pri- 
mary nuclei ratio. This component would have to be local in the sense of being specific to the solar neighbourhood, 
so that the well-known "Local Bubble" phenomenon is a natural candidate. 

The low-densi ty region around the Sun, filled with hot H I gas, is called the Local Bubble (hereafter LB) (e.g., 
Sfeir et allll999h . The size of the region is about 200 pc, and it is likely that it was produced in a series of supernova 



(SN) explosions, with the last SN explosion occuring approximate ly 1-2 Myr ago, or 3 SN occuring within the last 5 
Myr. Most probably its pr ogenitor was an OB star association (e.g., Ma iz-Apell anizl 2001 : Berghofer and BreitschwerdtL 



Myr. Most propaDly its progenitor was an UK star association (e.g., Maiz-Apellaniz, zUUt; cergnoter ana creitscnwer 
l2002h . lKnieetal.N l999h report about an excess of 60 Fe measured in a deep ocean core sample of ferromanganese 



crust sugges t ing a SN explosion about 5 Myr ago at 30 pc distance and the deposition of SN-produced iron on earth. 
ISonett et al.l (|1987) report an enhancement in the CR intensity dated about 40 kyr ago, which is interpreted as the 
passage across the solar system of the shock wave from a SN exploding about 0. 1 Myr ago. It could also be tha t 
"fresh" LB contributions from continuous acceleration in the form of shock waves ( Bykov and Fleishmarl 1992 ). 
and/or energetic particles coming directly from SN remnants still influence the spectra and abundances of local CR. 
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Fig. 1 . Left: Calculated p flux in a model with S = 0.47 and different normalization values D (xio 28 cm s 2 ). Solid 
curves - D = 3.3 at p = 3 GV, upper curve - local in terstellar (L I S), low er curve - modulated (3> = 550 MV). Dots 
- Dn = 2.6, dashe s - D = 4.3. Data: BESS 95-97 lOrito et all l2QQ0l) . BESS 98 lAsaoka et all l2002ft MASS91 
jBasini et all h 999ft . CAPRICE98 l Boez io et all |2001|> . Right: B/C ratio calculated with LB contribution, E b = 500 
MeV. The lines ar e coded as on the left panel . Lower curves - LIS, u pper - mo dulated (<& = 450 MV ). Data below 
200 MeV/n: ACE l Davi s et allli oOQ). Ulysses iDuVer nois et al.lll996bh. Vovaaer iLukasiak et "ail 1 1999): high energy 
data: HEAO-3 (Engelmann et al., 1990), for other references see Stephens and Streitmatter 11998). 



THE CALCULATION PROCEDURE 

We use the propagation model GALPROP 1 as described elsewhere ( Strong and Moskalenkol . 19981 : Moskalenko et al 



2002); for the present purpose the 2D cylindrically symmetrical option is sufficient. The spatial diffusion coefficient 
is taken as D xx = (5Dq{p/ po) s ; the diffusion in momentum space and other details can be found in our e arlier papers. 
The n ucleon injection spectrum of the Galactic component was taken as a modified power-law in rigidity Jjones et all 
200 lh . dq(p)/dp oc / \fl + (p/2)~ 2 , for the injected particle density. The heliospheric modulation is treated us- 



ing the force-field approximation. The LB spectrum is taken in the form dBvkov and Fleishmarl Il992l for continous 
acceleration by interstellar shocks): df /dp oc p~ v exp(— p/pb), where p is the rigidity, and p& is the cut off rigidity 
parameter. In terms of kinetic energy per nucleon E this can be re-written as df jdE = a(Z, A) A ^ E ^, m ^ p~ ri e~ p l pb , 
where a(Z, A) is the abundance of a nucleus (Z, A), Z, A are the nucleus charge and atomic number, correspondingly, 
m is the atomic mass unit, p is the momentum per nucleon, pb = 4 \/ (Ej, + m) 2 — rri 2 . We show the results obtained 
with rj = 1, but the particular spectral shape of the LB component is not important as long as it decreases sharply 
towards high energies and is much softer than the Galactic CR spectrum. 

The procedure to tune the CR elemental abundances, secondary/primary nuclei ratios, and p flux we adopted was 
as follows. The high energy part of B/C ratio and p flux measurements are used to restrict the parameters of the 
diffusion coefficient Do and 5, while the low energy part of B/C ratio is used to fix a value for the reacceleration level 
and define the parameters of the LB component. The Galactic CR elemental source abundances are tuned using the 
CR abundances at high energies where the heliospheric modulation is weak. The LB isotopic abundances are tuned 
to match the low-energy data by ACE and Ulysses. 

Figure n (left) illustrates the process of fixing the normalization of the diffusion coefficient using p's. The p flux 
shown is calculated with 5 = 0.47 and different normalizations in the diffusion coefficient, Dq = 2.6, 3.3, 4.3 x 10 28 
cm s -2 at p = 3 GV (for p's this corresponds to kinetic energy ~ 2 GeV). The injection index 7 is taken equal to 
2.28, and the Alfven speed va = 23 km s _1 . The p flux at maximum, ~ 2 GeV, appears to be quite sensitive to the 
value of the diffusion coefficient and allows us to fix it at Dq = (3.3 db 0.8) x 10 28 cm s -2 . The exact value of 5 is not 
critical since we compare with the p measurements at maximum. Inelastically scattered p's, the "tertiary" component, 
appears to be important at low energies onl y in the ISM. Figured (r i ght) shows corresponding calcul ations of the B/C 
ratio. A halo height of Zh = 4 kpc is used (Sjrong and Moskalenkol. 12001 : Moskalen ko et all 2001 ah . 

1 GALPROP model including software and data sets is available at http://www.gamma.mpe-garching.mpg.de/~aws/aws.html 
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The local interstel l ar (LI S) p and He spectra used in our p calculations are the best fit s to the data as des cribed 
i n iMoskalenko et ail J2002h . Because of the measurements with large statistics (BESS - Sanuki et al. 200(1 AMS 



Alcar az et aD boOO). d weak heliospheric modulation above 10 GeV, the error arising from uncertainties in the 
primary spectra is only ~ 5%. We note that as in the case of other nuclei, there should be an LB contribution to p 
and He spectra. The Galactic injection spectra of p's and He should thus be significantly flatter below several GeV to 
match the data points at low energies. This does not influence the p production because (i) the LB does not produce a 
significant amount of secondaries, and (ii) the p threshold production energy is high, ~ 10 GeV. 

Further tuning can be done using the high energy part of the B/C ratio, which is not influenced by heliospheric 
modulation and supposedly contains only a Galactic component of CR. Index 5 ~ 0.47 is chosen as giving the best 
match. A pure Kolmogorov spectrum, 5 = 1/3, thus seems to be excluded by p spectrum data taken in combination 
with B/C ratio data. The low energy B/C data were used to fix the LB spectral parameter at E f, = 500 MeV, but 
Eb = 400 — 600 MeV also provides good agreement with B/C data. In this way, we have been able to obtain a model 
simultaneously fitting p, He, p, and B/C data. 

ABUNDANCES IN COSMIC RAYS AND COSMIC RAY SOURCES 

The Galactic CR source elemental abundances are tuned ( at a nominal reference e nergy of 100 GeV), by a least 
squares procedure, to the abundances measured by HEAO-3 ( Engelm ann et al. I [l990h at 5.6. 7.5. 10.6. 16.2 GeV/n. 
Relative isotopic abundances at the source are taken equal to solar system abundances (Anders a nd Grevesse . 19891) . 
The key point in the fitting procedure is to obtain the correct abundance of B. 

The LB elemental abundances are tuned simultaneously 

with spectra using the low energy part of the B /C ratio 

and is otopic abundances at 200 MeV/n from ACE jw iedenb eck et 
2001) and at 185 MeV/n from Ulysses dDuVernois and Th aver. 
1996). For many elements ACE and Ulysses abundances 1 00 
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differ by 10%. For this reason, to the statistical errors we 
added 5% systematic error. The fitting procedure is influ- 
enced by the adopted value of the modulation potential. 
The following pairs of modulation potentials give almost 
the same final result: (ACE, Ulysses) = (400 MV, 700 
MV), (450 MV, 600 MV), (500 MV, 500 MV). The data 
deviate from calculations in both directions, which mean 
that we are unlikely to introduce essential systematic error 
by assuming a wrong value of the modulation potential. 

The diffusive reaccleration model with an LB compo- 
nent shows good overall agreement with data including 
secondary to primary ratios, spectra, and abundances. The 
derived Galactic CR source abundances and LB source 
abundances are pl otted in FigureElrelative to the solar sys- 
tem abundances I Grevesse a nd Sauval, 1998). The impor- 
tant result is that CR source and LB source abundances of 

all major elements (C, O, Ne, Mg, Si, S, Ar, Ca, Fe, Ni) are in good agreement with each other. Abundances of the Si 
group, Ca, Fe, Ni are near the solar system abundances. Abundances of other elements in Galactic CR and LB sources 
are mostly consistent with each other, and with solar system abundances, within a factor of 2. Relative to silicon, C, N, 
O, Ne, S are underabundant bot h in CR source and LB. This corresponds to the well-known first-ionization-potential 
or volatily correlation (see, e.g.. lMever et all 1 1998). 

Spectra of B, C, O, Fe are shown in Figure|3]for modulation levels 450 and 800 MV. In case of C, the normalization 
coefficient in the LB component is fixed as a(6, 12) = 6.35 x 10~ 4 cm~ 2 s _1 sr _1 for r\ = 1. Since the elemental 
abundances are tuned at both high and low energies, the calculated sub-Fe/Fe ratio also agrees well with data. 

Secondary nuclei K, Sc, Ti, V appear to be overabundant in the LB sources relative to the solar system (shown 
as upper limits in Figure 13, though the derived absolute LB abundances are not large. The derived LB abundance 
of Ti does not exceed that of Cr while the derived abundances of Sc and V are not larger than that of Mn. One 
possible reason for this excess is the uncertainty in the production cross sections, which is especially large for these 



Fig. 2: Derived abundance ratios Galactic-CR-source 
to solar-system (GCRS/SS) and LB-source to solar- 
system (LBS/SS), normalized to Si. Relative abun- 
dances for K, Sc, Ti, V are shown as upper limits. The 
dotted lines plotted at 1/2 and 2. 
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Fig. 3. Spectra of B, C, O, and Fe calculated with LB contribution. Upper curves - LIS, lower c urves - modulated 
using force field approximation ($ = 450 MV - solid curves, $ = 8 00 MV - dashes). Data: A CE j Davis et al., 2000, 
I2001|) . HEAO-3 (Enael mann et all Fl 990K for other references see Stephens and Streitmatter (1998). 



nuclei. Sometimes there is no measurement at all, while phenomenological systematics are frequently wrong by a 
factor of two or more. Often, there is only one measurement at ~ 600 MeV/n, which has to be extrapolated in 
both directions. This allows only a nearly flat Webber-type or Silberberg-Tsao-type extrapolations while the real 
cross secti ons usually have large reso nances below s everal hundred Me V, and decrease with energy above a few GeV 
(see, e.g., Moskalenko et al. , 2001 ah . We note that Davis et al. (2000) used semiempirical cross sections based on 



Webber et al. ( 1990) and also predicted fluxes of sub-Fe elements which are too low. 

An estimate of the overall error, which is reflected in the derived LB source abundances, can be obtained by 
assuming the absence of K, Sc, Ti, V in the LB source. Such a discrepancy is below 20%, and can be removed by 
allowing the production cross sections to increase at low energies by ~ 15 — 20%, which seems plausible. Another 
possibility is errors in flux measurements of the rare CR species. Ulysses and ACE measurements are not always in 
agreement. Note that even for such an abundant nucleus as Fe, which is the main contributor to the sub-Fe group, the 
discrepancy exceeds 10%, while the disagreement in abundance of Sc is ~ 30%. 

The derived source overabundance of sub-Fe elements in the LB could also in principle arise from composition 
differences between the ISM in the LB and solar or Galactic average ISM. This is suggested by the fact that the 
relative abundances of secondary elements in the LB sources are systematically larger than in the Galactic CR sources 
(Figure 

Isotopic distributions in cosmic rays 

Be and B isotopes are assumed all secondary, thus there is no possibility to tune them. Our calculation shows 
perfect agreement with the data on relative isotopic abundances of Be and B. This is in contrast with a "stan- 
dard" reacceleration model, wh ere we obtained a 15% discrepancy with relative abundances of Be and Be isotopes 
(iStrong and Moskalenkd. 12001). Abundances of stable isotopes of other elements are not very conclusive because 
they are present in the sources, but O and Si isotopic distributions still agree very well with data by ACE, Ulysses, 
and Voyager assuming only 16 and 28 Si isotopes are present in the LB component. C and N isotopic distributions do 
not agree too well, but this may point to a problem with cross sections. The calculated ratio 13 C/ 12 C ~ 0.11 at 120 
MeV/n ($ = 500 MV) in the model with LB con tribution is still a fac tor ~ 1.5 too large compared to the measured 



value , 0.0629 ± 0.0033 (Voyager 50-130 MeV/n, Webber et al., 1996) and 0.078 ± 0.011 (Ulysses 100-200 MeV/n, 



iDuVernois et al. . 1996ah . which may be connected in part with overproduction of 13 C by 15 N. If we replace the cu- 



mulative cross section 15 N+p — > 13 C with cross section 14 N+p ^ 13 C, the calculated ratio 13 C/ 12 C will be lowered by 
10% as estimated. Assuming the absence of the isotope 13 C in the Galactic CR sources gives another 10% reduction. 
Altogether these corrections yield 13 C/ 12 C ~ 0.09, close to the data. 

In case of radioactive isotopes Be, Al, CI, usually used as "clocks" in CR, the agreement with the most accurate low 
energy data by ACE is very good and all the ratios are consistent with each other indicating that a halo height z h = 4 
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kpc is a good estimate. Higher energy data by ISOMAX (Be) are also consistent with calculations considering the 
large error bars. The 54 Mn/Mn ratio indicates a somewhat smaller halo, but this may be related to uncertainty in its 
half-life and/or production cross section. The half-life of 54 Mn against 0~ decay is the most uncertain among the 
four radioactive isotopes - it is only one which is not measured directly. Apart from the half-life, a possible source of 
errors can be production cross sections of Mn isotopes . Only the reaction n at Fe+p — > 54 Mn on a natural sample of 
Fe has been measured well enough (see compilation bv lMashnik et al. . 1998h . 



CONCLUSION 

In a previous paper we have shown that the p flux and B/C ratio appear to be inconsistent with measurements when 
computed in standard diffusion/reaccelation models. In this paper we have demonstrated that this discrepancy can be 
resolved if some part of the local CR consists of a "fresh" component accelerated in the LB. The independent evidence 
for SN activity in the solar vicinity supports this idea. 

Combining the measurements of the p flux and B/C ratio to fix the diffusion coefficient, we have been able to 
construct a model consistent with measurements of important nuclei ratios in CR and derive elemental abundances 
in the LB. Calculated isotopic abundance distributions of Be and B are in perfect agreement with CR data. The 
abundances of three "radioactive clock" isotopes in CR 10 Be, 26 Al, 36 C1 are all consistent and indicate a halo size 
Zh ~ 4 kpc based on the ACE data. 54 Mn indicated a smaller halo, but this may be related to its half-life uncertainty 
and/or cross section errors. The derived fraction of LB component in CR is small compared to Galactic CR and has a 
steep spectrum with a cutoff above several hundred MeV/n. Other experimental data (except maybe the overabundance 
of Sc, Ti, V) do not contradict this hypothesis. The derived source overabundance of sub-Fe elements in the LB may 
be caused by trivial uncertainties in the production cross sections, or could in principle arise from composition and/or 
evolution differences between the ISM in the LB and solar or Galacti c average ISM. 



Recently there has appeared some indication (Huan g et all 12001) that the atmospheric contribution to the p flux 



measured in the upper atmosphere is underestimated by ~30%. This means that the flux of p's in CR in reality may 
be lower at the top of the atmosphere by at least 25-30%. If the latter is true, the reacceleration model (even without 
LB) could still be the best one to describe propagation of nucleon species in the Galaxy. 
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